The open-charm radiative and pionic decays of molecular charmonium 7(4274) 
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In this work, we investigate the decay widths and the line shapes of the open-charm radiative and pionic 
decays of ¥(4274) with the £),fl,o(2317) molecular charmonium assignment. Our calculation indicates that 
the decay widths of 7(4274) D*D'-y and K(4274) -> D+D-n° can reach up to 0.05 keV and 0.75 keV, 
respectively. In addition, the result of the line shape of the photon spectrum of y(4274) — > D^D'^y shows that 
there exists a very sharp peak near the large end point of photon energy. The line shape of the pion spectrum of 
K(4274) — » D'l.Dl^nP is similar to that of the pion spectrum of 7(4274) — » D'^D'^y, where we also find a wery 
sharp peak near the large end point of pion energy. According to our calculation, we suggest further experiments 
to carry out the search for the open-charm radiative and pionic decays of 7(4274). 

PACS numbers: 12.39.-x, 13.75.Lb, 13.20.Jf 



I. INTRODUCTION 

Recently the CDF Collaboration reported an explicit en- 
hancement structure with 3.1cr significance in the J/if/(f> in- 
variant mass spectrum of the — > K^J/i(/if> process. In 
this work, we refer to this new enhancement structure by the 
name 7(4274). Its mass and width are M - 4274.4+g^(stat) ± 
1.9(syst) MeV/c^ and T = 32.3!^[^(stat) + 7.6(syst) MeV/c^ 
ilt], respectively. We need to specify that this new result ap- 
pearing in the J/i^<p invariant mass spectrum is based on a 
sample of pp collision data at -y^ = 1-96 TeV with an in- 
tegrated luminosity of about 6.0 fb ' Ol. Additionally, the 
experiment presented in Ref. ^ also confirmed the observed 
y(4140) previously announced in Ref. 

Before finding the 7(4274) structure, there had been six 
charmonium-like states observed in B meson decays, which 
include X(3872) in B ^ J/i/zn^n-K id, 7(3940) in B ^ 
JHnjK iH H, F(4140) in -> J/if/(/) K ^ 



Z+(4430) in 

, and Z+(4015) and Z+(4248) in B Xci^^K 



where we use the underlines to mark the corresponding 
decay channels of charmonium-like states observed in exper- 
iments. The evidence of 7(4274) revealed by CDF ^ not 
only has made the spectroscopy of charmonium-like states ob- 
served in B meson decays abundant, but has also stimulated 
theorists' interest in revealing its underlying structure. Study- 
ing 7(4274) will improve our understanding of the essential 
mechanism resulting in these structures. Very recently, LHCb 
also confirmed the observation of 7(4274) with 53+ events in 
the B+ ^ JI\p<pK* decay H]. 

Since 7(4274) was observed in the Jl\l/<p invariant mass 
spectrum, we can conclude that the quantum numbers of 
7(4274) are = 0+(7++) with 7 = 0, 1, 2 if 7(4274) ^ 

JI4'4' occurs via S-wave. If explaining 7(4274) as a can- 
didate of charmonium, 7(4274) should be a P-wave state 



with the second radial excitation. In Ref. i^, the pre- 
dicted total widths of the second radial excitations of XcO 
and Xc\ are larger than the width of 7(4274). In addition, 
X(4350), which was observed in the yy — > (pJlif/ process 
with mass m - 4350.6+5j(stat) ± 0.7(syst) MeV and width 
r = 13^g-'(stat) + 4(syst) MeV fl^l, was explained as the can- 
didate of the second radial excitation of Xc2 1^, where the 
measured parameters of 7(4274) are different from those of 
X(4350). However, at present we cannot fully exclude the 
P-wave charmonium explanation of 7(4274), since the uncer- 
tainty of the quark pair creation model is not under control 

Besides the conventional charmonium assignment to 
7(4274), 7(4274) as molecular charmonium can be produced 
by B meson decay. The comparison of the existing ex- 
perimental information of charmonium-like states observed 
in B meson decays reflects a common property; i.e., these 
charmonium-like states are near the threshold of the corre- 
sponding charmed meson pair, which has provoked the in- 
vestigation of whether these observed charmonium-like states 
X(3872), Z+(4430), Z+(4015)/Z+(4248), 7(3930), 7(4140), 
and 7(4274) can be explained as the corresponding molecu- 
lar charmonia lfTTl - l2ll . The production of molecular charmo- 
nium via the B meson decay is allowed. The cc pair is created 
from the color-octet mechanism in the weak decays of the B 
meson; then c and c, respectively, capture q and q to form a 
charmed meson pair, where a color-octet qq pair is popped out 
by a gluon. Thus, a pair of the charmed mesons with low mo- 
mentum easily interact with each other to form the molecular 
charmonium ll20ll . 

As indicated in Refs. dHO], the mass of 7(4274) is near 
the threshold of DiZ)jo(2317), which is similar to the situa- 
tions of 7(4140) and 7(3930), since_ 7(4140) and 7(3930) are 
assigned as molecular states of Z)*D* and D*D*, respectively 
ll20l] . Thus, it is natural to deduce that 7(4274) can be an 
S-wave DjZ)jo(2317) molecular charmonium with the flavor 
wave function 
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17(4274)) = + 



(1) 
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which is also supported by the dynamical calculation of its 
mass when assuming F(4274) as an S-wave DjZ)io(2317) 
molecular state with spin-parity 7^ = 0" lfl6l I20I1 . 

Thus, studying 7(4274) under different assignments will be 
helpful in distinguishing its possible structure explanations. In 
this work, we mainly focus on the decay behavior of 7(4274) 
as a charmonium-like molecular structure. Performing the 
study of the decay behavior of 7(4274) can provide crucial 
information for testing this molecular assignment to 7(4274). 
Considering only the preceding reasons and the present theo- 
retical research status of 7(4274), in this work we investigate 
the open-charm radiative and pionic decays of 7(4274), which 
includes the calculation of the branching ratios and the study 
of the line shape of these decays. Our study will give hints for 
further experimental studies of 7(4274), especially in search- 
ing for other decay channels of 7(4274). 

This work is organized as follows. After this introduction, 
we illustrate the calculation details of the open-charm radia- 
tive and pionic decays of 7(4274) under the assignment of 
molecular charmonium. In Section |llll the numerical result 
will be presented. The paper ends with the discussion and 
conclusion. 



for the radiative and pionic open-charm decays of 7(4274) 
can be expressed as 

M[7(4274) ^ D+D'-y] 

= <Drrl'?^y|£>.o(2317)-)<£),o(2317)-D:i'Hr|7(4274)),(2) 
A4[7(4274) ^ Z)+D;7r°] 

= <Z);;r°|<H,|D,o(2317)-)<D,o(2317)-D:i'7/y|7(4274)> 
+<D>''|'K,|D,o(2317)+)<D,„(2317)+D;|'7/y|7(4274)), 

(3) 

where 'Ky describes the collapse of the S-wave DsDs[){23n) 
molecular state into and Dso(2317)". 'Hy or ']-(„ denotes 
the interaction of 0,0(2317)" with D^^y or D^n^. The main 
task of this work is to obtain <D,,o(2317)±Dj|'Hi|7>, which 
describes the collapse of the molecular charmonium 7(4274) , 
where we adopt the covariant spectator theory (CST) ||27| - |32|| 
to deduce it. 



A. Covariant Spectator Theory 



II. RADIATIVE AND PIONIC OPEN-CHARM DECAYS OF 

y(4274) 

Under the molecular charmonium assignment to 7(4274), 
it is interesting to investigate its radiative and pionic decays, 
since the photon and pion, respectively, from the radiative and 
pionic decays of 7(4274) can be easily detected in experi- 
ments. As a realistic research topic, both a theoretical esti- 
mate of their branching ratios and a study of the line shapes of 
the photon and pion spectra of the corresponding decays can 
reflect the internal structure of 7(4274) ll2iU2l . 

In this work, we focus on two groups of decay channels of 
7(4274), i.e., the radiative (D-^D*-y and D-D*^y) and pio- 
nic {D^D^tt'^) open-charm decays. In what follows, we take 
the decays 7(4274) ^ D^D*-y and 7(4274) -» D^D-ji^ 
as examples to illustrate the relevant calculation. Under the 
DjDjo(2317)" + /i.c. molecular state assignment, 7(4274)first 
dissociates into D, + 0^(2317)- or + Djo(2317)+. Then, 
the radiative decay 7(4274) — » D^D'fy and strong decay 
7(4274) -> D^D-n^ occur via the transitions D,o(2317)" 
D*-y and D,o(2317)- ^ D;7r°/D,o(2317)+ -> D+7r°, respec- 
tively, where D,o(2317)-/D,o(2317_}+ decays into D;7r"/'D>° 



by the 77 - n mixing mechanism 11241 - 12611 . The hadron-level 



descriptions of 7(4274) 
are shown in Fig.[T] 



D^D*-y and 7(4274) ^ D^D-jt" 





FIG. 1 : The radiative and pionic open-charm decays of 7(4274). 
The general expressions of the transition matrix elements 



The CST was proposed and developed to study the wave 
functions and the form factor of deuteron ||271 - |32| . where the 
CST is also referred to as the Gross equation. To some extent, 
the CST is an equivalent description of the Bethe-Salpeter 
(BS) equation when both are solved exactly. Thus, the CST is 
also widely applied to study two-body bound states. Since one 
particle is set to be on-shell, the Gross equation can be written 
easily in a form depending only on three momenta, which is 
different from the case of the BS equation. 

The Gross equation reads as 



(4) 



where 'V is the two-body interaction kernel and G is the prop- 
agator with particle 1 on-shell. In the CST, the Gross equation 
is shown in Fig. ID 




V 


+ 


V 










FIG. 2: (Color online.) Diagrammatic representation of the Gross 
equation for the two-body scattering matrix. Here "x" denotes the 
particle on-shell. 

Thus, bound state equations emerge automatically from the 
Gross equation. If the two-body system has a bound state at 
f 2 - one has 



|F)<F| 



(5) 
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where [F) is the vertex function and is finite at - M^. By 
Eqs. (HJi and (|5]l, we obtain 

ir> = vG|r>, (6) 

which corresponds to the description in Fig. |3] 

P/2 + P 

(1) 





P/'2 - P PI2- k 



FIG. 3: (Color online.) Diagrammatic representation of the Gross 
equation for the vertex function F. Here we use (1) and (2) to mark 
particle 1 and particle 2, respectively. 

After integrating over one obtains that the Gross equa- 
tion only depends on three-momentum 



r(p) 



J (27r)3 



^(p, k, W)G(k, W) F(fc), 



(7) 



where P - {W, 0) is the four-momentum of a two-body sys- 
tem, p and k are the relative momenta as shown in Fig. |3] 
Since particle 1 is on the mass shell, k - {ko, k), ko - 

Ei{k) - iW, and Ei{k) = ^m^ + k^. "Vip, k, W) is the inter- 
action kernel with particle 1 on the mass shell. 

For the case of 7(4274), the two-body propagator G{k, W) 
in Eq. ^ is 



G(fc, W) 



1 



4£i(fc)£2(fc) 
1 



1 



E2(k)- Ei{k) + W 



(8) 



E2(k) + Ei(k) - W 
The wave functions of the bound state can be introduced by 
1 



rip) 



np) 



V2£i(p)2£2(p)[£2(p) + Eiip) - W] 
1 



^J(2n)^2W 



np) 



V2£i(p)2£2(p) [Eiip) -Ei(p) + W] 



(9) 



(10) 



By the definition of wave functions and the equations for 
the vertex, we get the integral equations 

[E2ip) + Eiip)-W] ^\p) 

—r [Vip, k, W)i/,^ik) + V(p, k, W)rik)] ,(11) 

[E2(p)-Ei(p) + W] rip) 
<Pk 



[V{p, k, W)iA''(fc) + V(p, k, W)nk)] ,(12) 



where the potential is defined as 
V{p, k, W) 

1 



'V{p,k,W). (13) 



y/2Ei(p) 2E2(p) 2Ei(k) 2E2(k) 

The normalization of the wave function can be obtained by 
the normalization of the vertex. 



-! 



d 



1 = ( d'p r (P)^ [G(p, W)] Tip) - R 



(14) 



with 
R 



= J d'p j d'pr^(p)G(p,w)-^mp,p',w)] 



xG{p', W)T(p'). 



(15) 



Usually, the one-boson-exchange model (QBE) is widely 
applied to study hadronic molecular states |[l2l - [l8l 12011 . which 
makes the above calculation simple; i.e., the potential is not 
dependent on W, and = iH . 



B. Nonrelativized approximation of the CST 

In the following, we take the nonrelativized approximation 
in the CST, which was adopted in Refs [27-29] to study the 
wave function of the deuteron. Using the nonrelativized ap- 
proximation and the Fourier transform, the integral equations 
in Eqs. (fTTI) - (fT2l) can be transferred as 



-(y + e|i/'+(r) = -[V(r) + V(r)F(r)Vir)]il/^ir),{l6) 



(17) 



where F(r) - [m2-m\+W +V{r)Y'^ , Ei{p)+E2{p)-W ~ 
-e -H ^ with the reduced mass ^ and the binding energy of 
two-body system e-W-nii- 1112. 

Assuming mi 2 s> (V) as required from a loosely bound 
system, we obtain 



V- \ 

— + eU+(r) = -V{r)>]j*{r), 
M I 

rir) = 0, 



(18) 
(19) 



where Eq. ( fTSl l corresponds to the Schrodinger equation. The 
wave function in the momentum space can be obtained by the 
Fourier transform, i.e.. 



(2;r)3/2 /' 



(20) 



where p is the relative momentum. The wave function ij/^ip) 
satisfies 



d'p\r{pf^l, (21) 
from the normalization condition of the vertex in Eq. (Il4t . 
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C. Decay width 



III. NUMERICAL RESULT 



With the preceding preparation, we illustrate how to calcu- 
late the decay width of the radiative and pionic open-charm 
decays. The general expression of the decay width is 



dr 



1 

2E 



f 3 



|Mr(27r)V 



\ !=1 

2 



(2n)^2ei {Inflei {Infle^ 



(22) 



where q',(e,) (; = 1 , 2, 3) is the momentum (energy) of the final 
states. 

The decay amplitudes At for the sequential decay as shown 
in Fig. [T|can be written as 



M 



-^y(4274) -^D,,o(2317) 



'D^(2317) 



■q- 



(23) 



where !Ay(ati\) or ^d,o(2317) denotes the vertex for the de- 
cay of Y{4-274-) or Djo(2317). q and mD,„{23n) are the four- 
momentum and mass of the intermediate state D.5o(2317). By 
the formulism in Section HTAl and lirBl Eq. ( |23] ) can be further 
expressed in the center of mass frame of the decaying particle 
Fas 



y(4274) 



o(2317)(q) 



'Z5.^(2317) 



W - EoXq) + £'D,o(2317)(q) 

~ V2Mk(4274) ,/2(2;r)3 



mo. 



(24) 



where W is chosen as My{A2ia), which is the mass of 
F(4274). The above expression is a bridge connecting 
the decay amplitude and the wave function of the bound 
state. Thus, Eq. ( l24l i corresponds to the matrix element 
<D,o(2317)-D+|'Hy|y(4274)) in Eqs. Q-©. 

In addition, we also adopt the effective Lagrangians 

XD,„(2317)D:y = gy (8 ^.Aydf" fj^. - dyA^d^^l,.)(f>o,,o, (25) 
Xz5,„(2317)D,,;r = gjT<pD,<p7T<pDM (26) 

to depict the interactions of D5()(2317) with D*y and DjTt, 
respectively, where gy and g„ are the effective coupling 
constants. Thus, the amplitudes {D*"y|'7Yy|Dio(2317)") and 
<D;;r"|"K;,|D,o(2317)-> read as 

{D:-y\'Hy\D,o(23nr} = gy Ey,SD'Ag'"'qr ' qo- - ?;^.), 

(27) 
(28) 



(£)>0|<K,|Z),o(2317)^>=g,, 



where qyiqo',) and Syiso:.) are the momenta and the polariza- 
tion vectors of photon (£>j ), respectively. The results in Eqs. 
(|27]|-(|28T| correspond to ^Dm(23\i) in Eq. ( l23T l. 



In Ref. 112011 . the effective potential of 7(4274) was obtained 
by the OBE model 



V(r) = ""(r) = Vl^YiA,q'„m„r) 



with Y function 

y(A, K, m, r) = ■ 



Anr 



3 



1 



(29) 



(30) 



and ^1 - Vm- - k^, ^2 - VA^ - k^, q[^ - mo^oiim) - mo,, 
m,, = 548 MeV, and = 132 MeV, where 7(4274) is an S- 
wave 05^50(23 17) molecular charmonium. By the OBE po- 
tential, one finds the bound state solutions for 7(4274) shown 
in Fig. m where the dependence of the binding energy and 
the root-mean-square radius on different values of h and A are 
given (see Ref. |l20tl for more details). 
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FIG. 4: (Color online.) The obtained bound state solutions (binding 
energy E and root-mean-square r) of 7(4274) dependent on h values 
and cutoff A. 



For the decay of 7(4274) as shown in Fig. [T] in the cen- 
ter of mass frame of 7(4274), there exists a relation be- 
tween the four-momentum and three-momentum of the off- 
shell intermediate state D.5()(2317), i.e., q^ = ^y{4274) ~ 

2My(4274) ^TO^^ + q^, where ^/q^ also corresponds to the in- 
variant mass Mo^^y or Mo,7t- In Fig. |5] we list the variation 
of wave function ip^{q) obtained by the OBE model with 
when taking different h values, where we can reproduce the 
binding energy (My - mo^„{2}.ii) - ~ -11 MeV) with the 
corresponding A values. The results in Fig. |5] indicate that 
iff^(q) is not strongly dependent on li and A. 

Besides information on the wave function of the 
DsDso(23n) molecular state, the coupling constants gy and 
g„ can be extracted by the relations 



r(D;o(23i7)^Drr) = g 



^^^'"0^(2317) 



{qo'rqyfm 
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h=-1.6,A=1.35GeV 



h=-1.5,A=1.55GeV 




1.4 1.6 1.8 2 2.2 

V? (GeV) 

FIG. 5: (Color online.) The dependence of ifr^iq) on -s/q^. Here, -s/q^ 
corresponds to the invariant mass M^'y or M^^^ as shown in Fig.[Tl 



r{D%{23n)^Dy) = gl^-^ , (32) 

O„o(2317) 

where we use theoretical values to determine the constants, 
since these experimental partial decay widths are absent at 
present. In the literature Ii25, 26, 33-46J, the radiative and 
pionic decays of Z) ,0(23 17) were calculated under different 
structure assignments to Djo(2317). In this work, we take 
the typical values r(D,o(2317)- D'-jt) =10 keV iB and 
r(D,o(2317)- ^ D*-y) =1 keV iB, where D,o(2317) is as- 
sumed to be a charm-strange meson. We emphasize that our 
numerical result is dependent on the assumption on the par- 
tial decay width of D.so(2317), since different structure as- 
signments to Djo(2317) can result in different partial decay 
widths of Z)jo(2317). Discussing the structure of Z)jo(2317) 
is beyond the scope of this work. If the partial decay widths 
of Djo(2317) are determined in the future, the corresponding 
decay widths of 7(4274) open-charm radiative and pionic de- 
cays can be obtained by the present results multiplied by an 
extra factor. 

TABLE I: The radiative and pionic open-charm decay widths of 
¥(4274) with typical values of h and A. 



h 


A (GeV) 


Ty (keV) 


F;, (keV) 


-1.60 


1.35 


0.049 


0.74 


-1.50 


1.55 


0.050 


0.78 


-1.40 


1.85 


0.050 


0.75 


-1.30 


2.50 


0.049 


0.77 



In TableHl we give the radiative and pionic open-charm de- 
cay widths for 7(4274) with different combinations of h and 
A values. We find that these obtained partial decay widths of 
7(4274) are not sensitive to the corresponding combination of 
h and A. The calculated decay widths of 7(4274) D^D*-y 



and 7(4274) D'^D-tt'^ are around 0.05 keV and 0.75 keV, 
respectively. Of course, we admit that the above results are 
dependent on the decay widths of £)^()(2317) — > D* y and 
D^q(2317) — > £>^7r", since the coupling constants gy and g„ 
are determined by these theoretical values. We emphasize that 
in the current work the dissociation of 7(4274) is described by 
the wave function i^^(q). Thus, estimating the branching ra- 
tios of 7(4274) decays under the CST is a realistic approach, 
since the wave function adopted in the calculation is obtained 
by solving the Schrodinger equation with the OBE potential 
in Eq. (|29ll. 

In fact, we expect that the decays considered in this paper 
would be insensitive to details of the wave function. The rea- 
son is that the binding energy -11 MeV is much smaller than 
the range of forces, which is given by the mass of the eta in 
our model. This means that the size of the molecule is around 
1/(160 MeV). The decay channels through the decays of one 
constituent should "know" little about the short-distance be- 
havior of the wave function, since they are dominated by the 
long-distance part. This fact can be seen easily from Table 
HI where different sets of parameters give almost the same re- 
sults for the decay widths. One may further check Fig. |5] 
the long-distance part of the wave function is the same for 
different parameters, while one sees clear differences at short 
distances, which correspond to small values of q- and hence 
large values of q~. The observation of similar features shown 
in Fig. |6] should be strong evidence of an intermediate state 
with a well-defined mass of 2.3 GeV. 

Besides presenting the decay widths of 7(4274), we also 
carry out the study of the line shapes of the photon and pion 
spectra of 7(4274) D^D*--/ and 7(4274) D^D-jt^^ 
processes under the assignment of the DsDsQ{23n) molec- 
ular state to 7(4274). Here, we use the CERNLIB program 
FOWL to produce the Dalitz plots and the line shapes of 
the photon and pion spectra of 7(4274) — > D^^D*^y and 
7(4274) £)+D;7r" (see Fig. HJ. 

For 7(4274) D^D*j, an accumulation appears in the 
Dalitz plot with photon energy Ey ~ 0.19 GeV and Mo-y ~ 
2.3 GeV, which also reflects that D'*y and D^tt are from an 
intermediate state D.5o(2317). The line shape of the photon 
spectrum of 7(4274) — > D^D^y indicates that a very sharp 
peak exists near the large end point of photon energy, which 
directly corresponds to the accumulation in the Dalitz plot. 
The result of the line shape of the pion spectrum for the 
7(4274) — > D^D^TT^ process shows that a steep peak also ex- 
ists near the large end point of pion energy, which is similar 
to the situation of 7(4274) D^D*-y. 



IV. DISCUSSION AND CONCLUSION 

More and more observations of charmonium-like states X, 
Y, Z in B meson decays are providing us with a good plat- 
form to intensively study the internal mechanism for produc- 
ing these charmonium-like states, this is one of the most im- 
portant researc h top ics-full of challenges and opportunities-in 
charm physics 14711 . 

Stimulated by the recent evidence of 7(4274) Ulst] and the 
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FIG. 6: (Color online.) The Dalitz plot analysis and the photon and pion spectra for 7(4274) — > D^D* y and 7(4274) — > D^f , tt" decays. 
Here = -1.3 and A = 2.5 GeV. 



proposed DjZ)jo(2317) molecular charmonium explanation 
for y(4274) ll20ll . in this work we study the open-charm radia- 
tive and pionic decays 7(4274) -> D^^D* y and YiAllA) 
D^D^tt'^, where we not only present the calculation of its de- 
cays but also give the line shape of the photon and pion spec- 
tra of 7(4274) -> D^D*-y and 7(4274) -> D^D-jt" in detail. 
These theoretical predictions of the decay behavior of 7(4274) 
provide the information for further experimental studies of 
7(4274). 

We emphasize that 7(4274) — > D^D^n^ is peculiar 
to 7(4274) under the assignment of DjZ)jo(2317) molec- 
ular state, since this decay reflects the internal structure 
7(4274). 7(4274) D^D^^^^ occurs via the interaction 
of D,o(2317)" with D;7r" after the coflapse of 7(4274) into 
D+ and D,o(2317)-, where Djo(2317)" is a spe- 

cial decay channel observed by experiment ll24ll . Thus, the 
experimental search for a 7(4274) — > D^D^n'^ channel and 
the measurement of the line shape of the pion spectrum of 
7(4274) — > D^D'^n^ will be an interesting research topic. 

We must admit that the predicted partial decay widths of 
7(4274) are tiny compared with its total width, which makes 
the branching ratios of the considered decays about 10^^ ~ 
10"^, since the channels considered in this paper are either 



electromagnetic or isospin violation processes. Although 
these decay channels of 7(4274) cannot be the ideal ones ac- 
cording to the present status of CDF and LHCb, searching for 
open-charm radiative and pionic decays of 7(4274) could be 
considered topics for possible eff'orts in future experiments, 
such as the forthcoming Bellell IH] and SuperB ll49ll . 
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